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ABSTRACT 

The central parsec around the super-massive black hole in the Galactic Center 
hosts more than 100 young and massive stars. Outside the central cusp {R ~ 1") 
the majority of these O and Wolf-Rayet (WR) stars reside in a main clockwise 
system, plus a second, less prominent disk or streamer system at large angles with 
respect to the main system. Here we present the results from new observations 
of the Galactic Center with the AO-assisted near-infrared imager NACO and the 
integral field spectrograph SINFONI on the ESO/VLT. These include the detec- 
tion of 27 new reliably measured WR/0 stars in the central 12" and improved 
measurements of 63 previously detected stars, with proper motion uncertainties 
reduced by a factor of four compared to our earlier work. Based on the sample 
of 90 well measured WR/0 stars, we develop a detailed statistical analysis of 
their orbital properties and orientations. We show that half of the WR/0 stars 
are compatible with being members of a clockwise rotating system. The rotation 
axis of this system shows a strong transition from the inner to the outer regions 
as a function of the projected distance from Sgr A*. 



^Max-Planck- Institute for Extraterrestrial Physics, Garching, Germany 
''GRAAL-CNRS, Universit Montpelier II, Montpelier, France 
'^Department of Physics, University of California, Berkeley, USA 

'^Leiden University, Leiden Observatory and Lorentz Institute, , NL-2300 RA Leiden, the Netherlands 

''LESIA, Observatoire de Paris, CNRS, UPMC, Universit Paris Direrot, Meudon, France 

^Faculty of Physics, Weizmann Institute of Science, Rehovot 76100, Israel 

^Department of Physics & Astronomy, University of Leicester, Leicester, UK 

'^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, USA 

'School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 

*correspondence: H. Bartko, hbartko@mpe.mpg.de 



- 2 - 



The main clockwise system either is either a strongly warped single disk with a 
thickness of about 10°, or consists of a series of streamers with significant radial 
variation in their orbital planes. 11 out of 61 clockwise moving stars have an 
angular separation of more than 30° from the local angular momentum direction 
of the clockwise system. The mean eccentricity of the clockwise system is 0.36 ± 
0.06. The distribution of the counter-clockwise WR/0 star is not isotropic at the 
98% confidence level. It is compatible with a coherent structure such as stellar 
filaments, streams, small clusters or possibly a disk in a dissolving state: 10 out 
of 29 counter-clockwise moving WR/0 stars have an angular separation of more 
than 30° from the local angular momentum direction of the counter-clockwise 
system. The observed disk warp and the steep surface density distribution favor 
in situ star formation in gaseous accretion disks as the origin of the young massive 
stars. 

Subject headings: Galaxy: center - stars: early-type - stars: formation - stars: 
luminosity function, mass function - stellar dynamics 



Introduction 



The Galactic Center (GC) is a uniquely acces sible laboratory for studying the properties 



and evolution of ealactic nuclei ffor reviews see e.e. Genzel & Townes 1987 




Morris & Serabvn 


1996 




Mezeer et al. 


1996; 


Alexander 


2005). At a distance of about 8 kpc 


Eisenhauer et al. 


2003a 


Gillessen et al. 


2009; 


Ghez et al. 


2008; Groneweeen et al. 2008; 


Tripp e et al. 


2008 





processes in the Galactic Center can be studied at much higher resolution compared to any 
other galactic nucleus. Stellar orbits show that the gravitational potential to a scale of a few 
light hours is dominated by a concentrated mass of about 4 x IO^Mq. It is associated with 
the cor npact radio source Sgr A*, which mus t be a massive black h ole, beyond any reasonable 
doubt JSchodel et al.ll2002l : Ichez et al . I [2OO5I : I Gillessen et al.ll2009h . We will adopt a distance 
and a mass of Sgr A* of i?o = 8 kpc and M^gr a* = 4.0 x 10^ Mq for all analyses presented 
in this paper. The evolution and the star-formation history in the central pc of the Galaxy 
may also be used as a probe of star formatio n processes near supermassive b lack holes in 
general, also relevant to other galactic nuclei (jCoUin fc Zahnll2008l : iLeviru 120071 ). 



The central parsec of the Galaxy contains about a hun dred massive young stars. The 



majority are 0-type supergiants and Wolf-Rayet (WR) stars ( 


Forrest et al. 1987 


Allen et al. 


1990 


Krabbe et al. 


1991: 


Naiarro et a.. 


.|l994 


: Krabbe et al.lll995l: Blum et al. 


I995I; 


Tamblvn et al. 


1996 


Naiarro et al. 


1997 


Genzel et al. 


2003; 


Ghez et al.ll2003l: 


Paumarc 


et al. 


2006 


:lMartins et al. 


2007 


with an estimated aee of about 6 x 10^ vears. ICenzel et al. 


(199( 


I . 


2000. 2003): 
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Levin fc Beloborodovl (120031 ) ; iBeloborodov et al.l (120061 ) ; iPaumard et al.l (120061 ) inferred that 
most of the dynamical properties of the WR/0 stars (located at projected distances to 
SgrA* between 0.8" and 12") are compatible with belonging to either of two moderately 
thick counter-rotating stellar disks. However, iTanner et al.l (120061 ) were on ly able to assign 



15 ou t of 30 early-type stars near the Galactic Center as disk members. iLu et al.l (12006 



20081 ) confirm one stellar disk but do not observe a significant number of stars in the other 



one. 



The existence of these young massive stars indicates that star formation must have re- 
cently taken place at or near the Galactic Center within the last few million years. This is 
surprising, since regular star formation processes are likely to be suppressed by the tidal forces 
fro m the massiye black hole. Many scenarios h a ye been suggeste d for the origin of these stars 



see 



Alexander! (120051 ): iPaumard et al.l (120061 ): iPaumardl (120081 ) for recent reviews). These 



include in situ star formation through gravitational fragmentation of gas in disk(s) formed 
from infalling molecular cloud(s); transport of stars from far out by an infalling young stellar 
cluster, or through disruption of binary stars on highly elliptical orbits by the massive black 
hole; and rejuvenation of old stars due to stellar collisions and tidal stripping. The young 
stars observed in the inner R ~ 1" are less massive B-stars (so called the 'S-stars') and 
are lik ely to ori ginate f rom a different scenario t hen t he O and WR stars (e.g. 



Perets et al. 



( 120071 ). but see iLevinI (120071 ): iLockmann et al.l (120081 )). Here we discuss only our observa- 
tions of the O and WR stars outs ide the central 0.8" ( other populations of young stars in 



the GC are discussed elsewhere; iGillessen et al.ll2009t iMartins et al.ll2009l ). and interpret 



them in the context of the two leading formation scenarios, the infa l ling cluster (IGerhard 



2001 



McMillan fc Portegies Zwart 



Kim et al. 



2003 



2005 



2004 



Giirkan fc Rasio 



Genzel et al. 



2003; 



2005) and the in situ formation (iLevin fc Beloborodov 



2OO3I : IPorteries Zwart et al.l I2OO3I: iKim k MorrisI I2OO3I: 



GoodmanI l2003l : iMilosavljevic fc Loebl |2004J : iNayakshin fc Cuadra 



Paumard et al.ll2006h ~ scenarios. The infalling cluster and the in situ formation see- 



narios can be distinguished by diffe rent phase space distributions of the stars (see also 
Paumard et al.l l2006l : iLu et al.l 120081 ). Key observables are the number of disks, the frac- 
tions of disk and isotropic stars, the disk orientation, thickness, eccentricity and warp as 
well as the radial density of the stars and the stellar mass function. 

In the following we present the results of new observations of the Galactic Center with 
the adaptive optics (AO) assisted near-infrared imager NACO and the integral field spectro- 
graph SINFONI on the ESO/VLT. These include the detection of 27 new reliably measured 
WR/0 stars in the central 12" and improved measurements of previously detected stars, 
with proper motion uncertainties reduced by a factor of four compared to our earlier work. 
Based on a sample of 90 well measured WR/0 stars, we develop a detailed statistical analy- 
sis of their orbital properties and orientations. To this end, we use a Monte-Carlo technique 
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to simulate observations of a large number of isotropic stars with the same measurement 
uncertainties as present in the data. From these simulated measurements, we determine the 
probability of finding coherent dynamical structures against isotropic stars. We find strong 
evidence for the existence of a warped disk in the distribution of the clockwise rotating stars 
and a non-random structure among the counter-clockwise rotating stars, which is possibly 
an additional disk. We then analyze the properties of the stellar disks using both the 3D ve- 
locity information and the stellar positions. We discuss the implications of our observational 
results for models for the origin of the O and WR stellar population in the GC. 

This paper is structured as follows: First, we describe our observations, the data selec- 
tion criteria and present the properties of our data set in section [2j Thereafter, in section 
El we describe our simulations of the observations of isotropically distributed stars and disk 
stars. In section H] we introduce our analysis method to search for features in the star distri- 
bution. In section 5 we present our results, including a thorough study of the significance of 
the counter-clockwise system, the determination of the orbital properties of the disk stars and 
a comparison to previous work. After a discussion of our results in section 6 we summarize 
our conclusions in section 7. 



2. Data 



2.1. Observations 



The data set previously analyzed by IPaumard et al.l (120061 ) contained 63 reliably (la- 
beled "quality 2") measured WR/0 stars in the innermost 12" and several candidates. In 
2006-2008 we carried out new observation s with the integral field spectrograph SINFONI 
( lEisenhauer et al.ll2003l : iBonnet et al.l 120041 ) at the ESO/VLT. We covered two regions west 
and north of Sgr A* with the AO scale (25 mas/pixel) resulting in a final K-band full 
width at half maximum (FWHM) of typically about 100 mas. In addition, we observed 
sixteen 4.2" x 4.2" fields with the 100 mas/pixel scale resulting in typical K- band FWHMs 



of about 200 mas. For some of th e fields we used the laser guide star facility (IRabien et al. 



20031 : iBonaccini Calia et al.ll2006l ). The location of the observed fields is indicated by black 
squares in figure [H T he details of the observations and the data analysis will be presented by 
Martins et al.l (120091 ). These observations resulted in the reliable detection of 27 new WR/0 
stars near the Galactic Center. 25 out of the 27 new sta rs are located at projected distances 
between 5" and 12" . Four of the new stars were listed by lPaumard et al.l (120061 ) as early type 
candidates (quality and 1). We determine d the radial velocitie s of these new stars by fitting 
the observed spect ra with template spect ra (IMartins et al.ll2007l ). We also updated the radial 
velocities given by lPaumard et al.l (120061 ) for all stars in the re-observed SINFONI fields. In 
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Fig. 1. — The sample of 90 WR/0 stars (niK < 14 and A.{vz) < 100 km/s) in the central 0.5 pc 
of our Galaxy: Blue circles indicate clockwise orbits (61 WR/0 stars) and red circles indicate 
counter-clockwise orbits (29 WR/0 stars). The black circles show projected distances of 0.8", 3.5", 
7" and 12" from Sgr A *. Squares indicate the exposed fields with SINFONI in the 25 mas/pixel and 
100 mas/pixel scale. The whol e inner 0.5 pc region is contained in lower resolution (250 mas/pixel 
scale) SINFONI observations iPaumard et al\\200a) . 



addition, for all these early-type stars we derived proper motions from the NAOS/CONICA 
( jRousset et al.ll2003l : iHartung et al.ll2003l ) imaging data set of the Gala ctic Center coveriii g 
six epochs from May 2002 to March 2007 in the 27 mas/pixel scale ( iTrippe et al.ll2008l ). 
Table [U summarizes the K magnitudes, positions, proper motions and radial velocities of the 
90 WR/0 stars used in our analysis. 



2.2. Data Selection Criteria 



In this work we focus on the analysis of the dynamics of WR/ O stars and B supergiants. 
The brig htest of the so-called S-stars within th e central 0.8", S2, is an early B dwarf (BO - 
B2.5 V) jMartins et al.lboOSah with tuk = 14.0 JPaumard et al.ll2006[ ). We define thk = 14.0 
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Table 1: K magnitudes, positions, proper motions and radial velocities of the 90 WR/0 stars 
used in our analysis. The full machine-readable table is available online. 



Star 


mK^ x(")^ 






o-(vra) 


vdec 


O'(vdEc) Vr cr(Vr) 


1 


13.7 -0.776 


-0.2771 0.822 


83.2 


0.7 


-57.1 


0.9 -75 26 



"The uncertainty on the observed magnitude is 0.1 mag. 

''Positions are relative to SgrA*, typical position uncertainties are as low as 0.0002". 
■^AU velocities are in km/s assuming Rq — 8 kpc. 



as the border between O and B dwarfs and only include early-type stars with uik < 14 
in our analysis. Hence we exclude B dwarfs — S-stars type stars — from the present 
analysis. 0/WR and B supergiants all have initial masses larger than ~ 15M0, while B 
dwarfs are less massive. The S-stars have different dynamical properties than the disk stars . 



They have isotropic orbits with large eccentricities (lEisenhauer et al.ll2005l : iGhez et al.ll2005 



Gillessen et al.ll2009l ). Our analysis shows that a large fraction of the B dwarfs in the region 
of the disks have different kinematics than the WR/0 stars. A thorough analysis of the 
properties of B dwar fs in the region of the WR/ O stars [R > 0.8") will be presented elsewhere 
jMartins et al.ll2009h . 



In order to perform a reliable analysis of the stellar dynamics, we require the observation 
of a high quality early-type spectrum, such that the error of the radial velocity is 100 km/s 
or smaller; A{vz) < 100 km/s. 

After applying these cuts our sample comprises 90 WR/0 stars. There are about 10 
candidate WR/0 stars with either a too large radial velocity uncertainty or for which no 
reliable proper motions could be determined due to crowding. No WR/0 star is reliably 
measured in our data set in the central 0.8" and beyond 12". Figure [T] shows the locations of 
these stars relative to Sgr A*, and indicates whether these stars are on clockwise or counter- 
clockwise orbits. We only covered a relatively small area beyond 12" by deep integral-field 
spectroscopic observations as indic ated in figure [H The la rge "frame" between 15" and 20" 
around Sgr A* shown in figure 1 of iPaumard et al.l (120061 ) did not contain any WR/0 star, 
but also had a poor signal to noise ratio. 
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Fig. 2. — The sample of 90 WR/0 stars in the central 0.5 pc of our Galaxy: (Left) Distribution of 
the statistical velocity errors. The distributions of velocity uncertainties in the x and y directions 
have both a mean of 5 km/s and an RMS of 3 km/s. The distribution of radial velocity uncertainties 
has a mean of 51 km/s and an RMS of 25 km/s. (Right) Distribution of K band magnitudes. 



2.3. Properties of the Data Set 



Figure [2] (left panel) presents the distribution of the statistical velocity uncertainties 
in the x, y and z directions. The distributions of velocity uncertainties in the x and y 
directions both have a mean of 5 km/s and an RMS of 3 km/s. Th e systematic uncertainty 
of the astrometric reference frame is 6.4 km/s (ITrippe et al.ll2008l ). The uncertainty of the 
distance to Sgr A*, which we used to convert the measured proper motions on the sky 
to velocities in km/s, introduces a systematic error of about 6% in and Vy. The mass of 



Sgr A* and its distanc e are tightly correlated: Msgr a* = (3.95±0.06) x 10''(i?o/8 kpc)^-^^MQ 
( IGillessen et al.ll2009l ). In order to evaluate the systematic errors introduced by the choise 
of i?o = 8 kpc and Msgr a* = 4 x lO^M© we ran all analysis steps also with Rq = 7.5 kpc, 
Msgr A* = 3.5 X IO^Mq and Rq = 8.5 kpc, Mggr a* = 4.5 x lO^M©. 



The combined statistical and systematic proper motion uncertainties in the lPaumard et al. 



(120061 ) analysis had a mean of 35 km/s and an RMS of 9 km/s. The reduced errors of the 
present proper motion data are due to the larger data set, a correction of the geometric 
image d istortion and smal ler systematic uncertainties of the movement of the coordinate 
system (ITrippe et al.ll2008l ). The distribution of radial veloc ity uncertaintie s has a mean of 
51 km/s and an RMS of 25 km/s, about the same as for the iPaumard et al.l (120061 ) data set 
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Fig. 3. — The sample of 90 WR/0 stars in the central 0.5 pc of our Galaxy: (Left) Distribution of 
projected distances to Sgr A* for all WR/0 stars (black full histogram), clockwise moving WR/0 
stars (blue dashed histogram) and counter-clockwise moving WR/0 stars (red dotted histogram). 
(Right) projected and normalized angular momentum on the sky j = Jz/Jz, max o,s a function of 
projected distance to Sgr A *. The error bars for most of the stars are smaller than the markers. 
Clockwise orbits correspond to j > and counter-clockwise orbits correspond to j < 0. 



(mean 46 km/s and RMS 26 km/s) but for 90 instead of 63 high quahty stars. Figure [2] 
(right panel) shows the distribution of observed K band magnitudes. 

For some stars with low projected distan ces from Sgr A*, significant accelerations are 



measured, which allowed iGillessen et al.l (120091 ) to determine full orbital solutions (see below 
figure [1] 

However, a timeline of 



Trippe et al. 


2008: 


Lu et al. 


2008) 



value of the z-coordinate (see 
5 years of NACO data has not yet been sufficient to constrain the 2;-coordinate for more 
than a handful of stars. We have chosen not to include these limits as priors in our analysis 
in order to have a data set with uniform errors independent of projected distance to Sgr A*. 
Instead we compare below the statistical results of the entire data set with the properties of 
the subset of stars with measured orbits. 

Figure [3] (left panel) shows the distribution of projected distances to Sgr A* for the data 
set of 90 WR/0 stars. The sample of WR/0 stars is not complete, i.e. the probability to 
detect a star with a given magnitude in an image (photometric completeness) and to identify 
its spectrum as early-type (spectroscopic completeness) is below one. The completeness 
depends on the respective field as well as on the apparent magnitude of the star (see e.g. 
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Table 2: Numbers and fractions of spectral subclasses of the 90 selected WR/0 stars. 





clockwise 


counter-clockwise 


Type 


Number 


fraction 


number 


fraction 


OB 


42 


0.69 


17 


0.58 


Ofpe/WN9 


5 


0.08 


4 


0.14 


WN 


8 


0.13 


2 


0.07 


wc 


6 


0.10 


6 


0.21 


sum 


61 


1.00 


29 


1.00 



Martins et al.l 120091 ). O stars in the magnitude interval = 13 — 14 typically have a 
combined photometric and spectroscopic completeness of 75%. To obtain the most reliable 
radial distribution a correction for completeness will eventually have to be applied. However, 
considering the emphasis of the current paper on the angular momentum distribution this 
correction is unnecessary. Figure [3] (right panel) shows, for each of the e arly-type stars, the 
projected and norma lized angular momentum on the sky j = Jz/ Jz,ma.x (IGenzel et al.ll2003l : 



Paumard et al.ll2006l ) as a function of projected distance to Sgr A*. There are 61 WR/0 stars 
on clockwise orbits (j > 0) and 29 WR/0 stars on counter-clockwise orbits (j < 0). Most 
of the WR/0 stars with projected distances below 3" are on projected clockwise tangential 
orbits (j ~ 1). For larger projected distances there are two concentrations of stars with j — 1 
and j ~ —1 (project ed clockwise and couri ter-clockwise) tangential orbits. The "diagonal 
feature" observed by iPaumard et al.l (120061 ) looks less pronounced in our data. Our data 
rather suggest two systems of stars with j — 1 and j ~ — 1 and a background of random 
stars. 

Table [2] lists the numbers and relative fractions of the different sub- types of the selected 
90 WR/0 stars. The fraction of different spectral sub-types is strikingly similar for the 
clockwise and the counter-clockwise rotating stars. This resemblance in content of massive 
stars strongly suggests that the clock wise and counter-clockwise stars are, with in measure- 
ment accuracies of ~ 2 Myr, coeval (IGenzel et al.l l2003l : iPaumard et al.ll2006l ). Moreover, 
there is no significant difference in the fraction of spectral sub-types with distance to Sgr A*. 



3. MC Simulation of Signal and Background 



To assess the probability of the observed features in the stellar distribution being com- 
patible with an isotropic star distribution, we simulated measurements of isotropically dis- 
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tributed stars and disks. For these simulations we assumed bound orbits, which are described 
by the following orbital elements (see appendix |X] for a visual representation of the chosen 
coordinate system and the orbital elements): 



1. Longitude of the ascending node, inclination 

■vv- 9j, (f)j\ Orientation of the orbital plane / the orbital angular momentum vector 

2. a: semi-major axis 

3. e: eccentricity 

4. uj: argument of periapsis 

5. r: time difference between periastron and the current position of the star. 



To simulate an isotropic cusp stellar distribution we used the following algorithm: 



1. Generate the direction of the angular momentum vector J /\J\ uniformly distributed 
on a sphere. Calculate VL and i. 

2. Generate uj uniformly in the interval [0, 27r] and r uniformly in the interval [0,T]. 

3. Generate e according to the distribution dA^stars /de oc e (Binnev &: Tremainj 1987 ) for 
isotropic orbits. Note that the clockwise disk of iPaumard et al. J2006h has a different 
distribution of eccentricities. 



4. Generate a according to the probability density dA^^stars/da oc a^^^"^ with /3 = 2 in 
the interval [0.2" < a < 40"]. We motivate this choice as follows: For a power-law 
cusp the radial star nur nber density is related to the surface nu mber density S(rdisk) = 



dAi'stars/dAdisk oc T (ISchodel et al.l l2003l : lAlexandeii l2005l ). raisk denotes the 3D 
distance in the plane of the disk. Each projection (e.g going from 3 D density to 2D 



(surface) density) decreases the power of the radial coordinate by 1. IPaumard et al. 



( I2OO6I ) measured the surface number density of the early-type stars in the disks as 



S(rdisk) oc and in our work we find S(rdisk) oc ?"dkk^^^°'^^5 see section 15.1.11 



For the simulation of a thick disk we assumed the following distribution of the stellar 
angular momentum directions: 



dS 



oc exp 
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S is the solid angle in which the stellar angular momentum points, tjj is the angular distance 
between the generated angular momentum direction of the star {(pj, 6j) and the disk angular 
momentum direction ((/>disk; ^disk)- In the following we will call the parameter the two- 
dimensional Gaussian sigma thickness. 

The generation of e and a in steps 3 and 4 implies a strong prior to the simulated 
isotropic stars drawn from a /3 = 2 cusp distribution. However, we will only use the simulated 
isotropic stars to determine expectation values of the mean and RMS of the distribution of 
stellar angular momenta per solid angle. To first order, these distributions depend only on 
the direction of the orbital plane (given by fl and i). In the simulations we only considered 
the potential of S gr A*. For simplici ty, we neglect in our simulations the mass of the cluster 



of late- type stars (iTrippe et al.ll2008l ). about 2 x IO^Mq enclosed in the innermost 10". This 
is justified as we are only interested in the distribution of the stellar angular momenta of the 
simulated stars, which is, to lowest order, independent of the enclosed mass. 

For each simulated star we calculated the true positions and velocities x, y, z and 
Vx-,Vy,Vz from the generated orbital elements. We simulated the measurement uncertainties 
for the position and velocities by adding random numbers to these positions and velocities 
following the distribution of errors in the measured data, see figure [2l 



4. Analysis Method to Search for Features in the Star Distribution 

The prime criterion for distinguishing a disk from an isotropic stellar distribution, which 
we assume as the null hypothesis, is the presence of a common direction of the angular mo- 
mentum vectors for all stars. The angular momentum vectors of isotropic stars are uniformly 
distributed over a sphere. The analysis requires three steps: 



1. computation of a density map of angular momentum vectors of the n observed stars 

2. computation of the mean and RMS density maps of angular momentum vectors for n 
MC simulated isotropic stars. 

3. computation of a significance map by comparing the density map of step 1 to the mean 
and RMS expectations derived in step 2. 
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4.1. Computation of a Density Map of Angular Momentum Vectors 

The computation of a density map of angular momentum vectors for sets of measured 
or MC generated stars proceeds in three steps: 

1. MC generation of 1000 values for the unknown 2;-coordinate for each star, computation 
of the angular momentum direction 

2. computation of the total reconstructed angular momentum distribution on a sphere 

3. computation of the density of reconstructed angular momentum vectors in a fixed 
aperture. 



4- 1.1. Generation of z- Coordinates 



The orbit of a star in a given gravitational potential is fully described by six orbital 
elements. For most stars there are only five measurements: the projected positions on the 
sky X and y, the proper motion velocities Vx and Vy, as well as the radial velocity Vz 0. Under 
these conditions, the direction of the angular momentum of any given star cannot be uniquely 
determi ned. All possible directi ons of the angular momentum lie on a one- dimensional curve 
(see e.g. lEisenhauer et al.ll2005l ). Assuming that the star is bound, one can derive an upper 
limit to the line-of-sight position of the star, which in turn limits the curve of possible angular 
momentum directions. 

The specific angular momentum (angular momentum per mass) J is given by: 

J = fx V . (2) 

As the 2;-coordinate is unknown, the angular momentum can only be determined as a function 
of unknown z: J = J{z). For a cluster in equilibrium the distribution of the 2;-coordinates is 
related to the star surf a ce nu r nber density dist r ibutio n E(rdisk) = dAgtars/d^disk, see section 
15.1.11 and iGenzel et al.l (120031 ) ; iPaumard et al.l (120061 ) . Assuming a power law of the surface 
number density dAgtars/d^disk oc (x^ + y'^)~^^'^, this density distribution tran slates into a 
density distribution for the 2;-positions of the measured stars (lAlexanderl |2005| ) : 



dA«. 



dz 



dA. 



oc 



dV 



xo,yo 



oc [X, 



+ 1/0 + ^') 



-(/3+l)/2 



(3) 



where Xo,yo are the measured positions on the sky. 



For the definition of the coordinate system used, see appendix \K\ 
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Fig. 4. — Input (black points) and reconstructed distributions of the semi-major axis and the time 
difference between periastron and the current pos ition of the sta r. The uniform z prior is shown by 
a red histogram, the uniform acceleration prior \Lu et al\\200di) is denoted by the green histogram 
and the z-prior used in this work by the blue histogram. 



For each star we generate 1000 values, assuming P = 2, according to section [5.1.11 and 



Paumard et al.l (120061 ). The choice of (3 is again a strong prior to our analysis. However, the 
generation of the z-values is based on the same stellar surface density distribution as the 
generation of the MC simulated stars (see section [3]). We also ran our analysis simulations 
with /3-values of 1.5 and 2.5 and obtained similar results. The limit |2;|max > \z\ is given by 
the requirement that the star must be bound to Sgr A*. 

The way we determine the distances along the line-of-sight is very important as concerns 
possible biases. 

As a consistency check of our generation method for z we compare the reconstructed 
orbital element distributions for MC generated isotropic stars with the input distribution into 
the MC simulation. Moreover, we have computed the dist ributions of orb ital elements for the 



case of uniform z priors and uniform acceleration priors (iLu et al.ll2008l ). The distributions 
of eccentricity and the argument of the periapsis match very well the input distribution for 
all three studied 2;-priors. Figure H] shows the input and reconstructed distributions of the 
semi-major axis and the time difference between periastron and the current position of the 
star. The uniform 2;-prior gives strong biases to reconstructed positions near the pericenter 
of the star. Also the uniform acceleration prior biases the reconstructed star position to the 
pericenter. Although the reconstructed distribution with our adopted z-prior shows small 
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Fig. 5. — cylindrical equal area projection of the distribution of the direction of the angular mo- 
mentum vector (cosj vs. for the star IRS16 CC for 1000 generated z-values. The z-coordinate 
is generated according to equation ^ ( assuming a power law stellar surface density with exponent 
(5 = 2). Only z-values of hound orbits are considered. Clockwise orbits have angular momentum di- 
rections in the upper hemisphere (cos 9 > 0) and counter-clockwise orbits have angular momentum 
directions in the lower hemisphere (cos 9 <0). 



biases, these are smaller than the biases for the other two studied z-priors. 



4.1.2. Angular Momentum Distribution on a Sphere 



Having calculated the angular momentum for the 1000 generated 2;-values, we can com- 
pute the distribution of directions of the angular momentum vectors J{z) on a sphere. We 
parameterize the sphere with the usual spheri cal coordinates (t> and 6. For our further analy- 



sis, we use a cylindrical equal area projection (lGalllll885l : |Peterslll983l ) of this sphere: (cos 6*7 
vs. We divide the projected map into a grid of 90 x 90 bins in and cos^. As an 

example, figure [5] shows a cylindrical equal area projection of the distribution of angular 
momentum vectors directions J{z) {cosOj vs. (pj) for the star IRS16 CC. The shape of the 
distribution of reconstructed angular momenta is independent of Msgr a*- It determines the 
maximum |z|-value for bound orbits and thus the minimum | cos^j|. 



- 15 - 



4-1.3. Computation of Angular Momentum Density Maps 

In a last step, the average density pa{0,(j)) of the reconstructed angular momentum 
directions is calculated over an aperture of a = 15° radius centered on the sky direction 



1000 I 1 /l[J{z),{6,(j))] < a 

p.{e, 0) = , . ^ ■ (4) 

1000 • 27r(l — cos a) 

This averaging maximizes the signal (mean integrated angular momenta over the aper- 
ture for a disk) to noise (RMS integrated angular momenta over the aperture for isotropic 
stars) ratio for a moderately thick disk with a two-dimensional Gaussian sigma of 10°. We 
chose a simple aperture in order to have well defined sky regions with a flat weight. A 
matched filter may yield even higher signal to noise ratios. Bins of the angular momentum 
density sky map, which are separated by less than the aperture size, are correlated. 



4.2. Density Maps of Angular Momentum Vectors for Simulated Stars 

Our null hypothesis is that the WR/0 stars are isotropic. In order to be able to 
reject this null hypothesis we have to accurately characterize both the average and the 
fluctuations of the distribution of isotropic stars. In order to show that the observed stars 
are compatible with a disk of stars we only need to characterize the average distribution of 
disk stars. Therefore, we simulated a total of 4 x 10^ isotropic stars and several disks of 
4 X 10^ stars each with different orientations and thicknesses. As an example, we selected 
isotropic stars with projected angular distances from Sgr A* between 0.8" and 12". For 
these stars we computed sky distributions of reconstructed angular momentum densities for 
a fixed aperture of 15° radius. Figure [6] (left panel) shows a cylindrical equal area projection 
of the average reconstructed angular momentum density. The right panel of figure E] shows 
the same distribution but for a dis k with a two-dimension al Gaussian sigma thickness of 10° 



oriented like the clockwise disk in iPaumard et al.l (120061 ) in the same radial interval (right 



panel). In case of isotropic stars, the sky distribution of reconstructed angular momenta is 
rather flat (e.g., ±13% variation). There is a small depletion close to cos 6* = caused by the 
choice of simulated stars with projected distances in the interval 0.8"-12". This projection 
along a cylinder causes a small bias of the actual ^-coordinate of the star to have a larger 
absolute value than the generated ^-coordinates. In our example the distribution of MC 
generated isotropic stars has a mean of zero and an RMS of 5.7", while the distribution of 
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Fig. 6. — cylindrical equal area projections of the sky distributions of the average density (fixed 
aperture of 15° radius) per star of the reconstructed angular momentum directions for isotropic 
stars with projected angular distances from the Galactic Center between 0.8" and 12" (left panel) 
and for a stellar disk oriented like the clockwise disk in \Paumard et al\ \200a ) in the same radial 
bin (right panel). The sky distribution of the average density of reconstructed angular momenta is 
rather flat in the case of isotropic stars. For a disk there is a peak at the simulated disk angular 
momentum direction. The simulated disk position is marked with the black circle. The plots are 
based on a simulated data set of 4: x 10^ isotropic stars and 4 x 10'* disk stars, see text. 



the reconstructed z values has an RMS of 5.4". This causes the small bias towards face-on 
reconstructed orbits. For a disk there is a peak at the simulated disk angular momentum 
direction. 

In the case of 90 isotropic stars the MC simulations predict an average v alue of 1.6 stars 



per 15 ° radius aperture centered on the direction of the clockwise disk in iPaumard et al. 



(120061 ). This value is similar to the zeroth order expectation of 90 x 47r/ (27r(l — cos 15°)) = 
1.53 stars. The expected RMS is 0.55 stars per 15° radius aperture. A simulated thick disk 
of 90 stars rather yields a reconstructed angular momentum density of 17.2 stars per 15° 
radius aperture. Hence a disk can be well distinguished from isotropic stars. 



4.3. Computation of Significance Maps 



The goal is to test whether the features of a given population of n stars are statistical 
fluctuations of an isotropic stellar distribution or not. We applied the same cuts in projected 
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Fig. 7. — Distribution of reconstruc ted angular momentum directions per 15° aperture centered at 
the position of the clockwise disk o nPaumard et al\ \200a ) for sets of 90 simulated isotropic stars 
(black histogram). The blue line shows a Gaussian fit. 



distance to the MC simulated stars as to the data. We grouped the simulated isotropic 
stars in N sets of n stars. For each set of simulated stars we computed a sky map of 
the density of reconstructed angular momentum directions per 15° aperture, again in a 
cylindrical equal area projection cos^ vs. with a 90 x 90 grid of bins. We histogrammed 
the distribution of N values for the density of reconstructed angular momentum directions 
per 15° aperture for each of the 90 x 90 bins. These distributions can be well approximated 
by Gaussian distributions. We produced two sky maps containing the mean and the RMS 
of these distributions. As an example, figure [7] shows the distribution of reconstructed 
an gular momentum direc tions per 15° aperture centered at the position of the clockwise disk 
of IPaumard et al.l (120061 ) for sets of 90 simulated isotropic stars. This distribution can be 



approximated by a Gaussian. 

We computed significance maps from the sky map of the density of reconstructed angular 
momentum directions of the ob served stars. W e define the significance for each bin of the 



sky map as (see equation 10a of iLi fc Mai (1l983l )): 



observed density - mean density of isotropic stars 

signincance = ^^^^ ^ — . (5) 

RMS density of isotropic stars 

As the expected distribution of angular momentum densities for isotropic stars is approxi- 
mately Gaussian, we will call an excess of x times the RMS density over the mean density 
expected for isotropic stars an excess of x a. 
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Fig. 8. — cylindrical equal area projection of the probability that the star IRS16 CC is part of a 
thin disk with an angular momentum direction given by {6,(j)), without assuming any prior on the 
z -position of the star. 



4.4. Definition of for Disks 

With the method described above we can determine the position and significance of 
possible disk features in the distribution of the early-type stars. The next step is to quantify 
if a given star is a candidate member of the disk. We define a ^^(^disk, 0disk) value that a 
star with measured 3D velocity (f^.m, "^^j/.m, "i^z.m) and 2D position {x^,ym) has an angular 
momentum direction n= J/\J\ = (cos0disk sin 6'disk, sin 0disk sin 6'disk, cos6'disk) by minimizing 
the following function with respect to (a, e, u, r): 



X^(6'disk,0disk) = Min 



+ 



Xm-a:(a,e,cij,r,Qdisk,»disk) 



+ 



+ 



Vy,in-Vyia,t,U>,T,Cldisk,idiak) 



Avy 



-v^{a,e,uj,T,ndiBk,idiBk) 
At;,.,n 



(6) 



This definition does not use any priors except bound orbits to Sgr A*. The orbital el- 
ements (a, e, tu, r, f2disk, "^disk) define the motion of the star (see section [3]). ^^disk and idisk 
are fixed by the disk angular momentum direction: i^disk = arctan(— tan0disk), "^disk = 
arccos(— sin 0disk sin ^^disk)- The above defined has only observational errors in the denom- 
inators, and thus it has no bias to any models or priors. Contrary to the definition of 
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Fig. 9. — Distribution of the reconstructed angular deviation ■0 of the disk stars from the common 
disk angular momentum direction for four MC generated disk of 1000 stars each, with 0° , 5°, 10° 
and 15° thickness. 



Levin fc Beloborodovl (120031 ) our definition takes both the velocity and the projected position 
of the star into account. As an example, figure [H] shows a cylindrical equal area projection 
of the probability (one degree of freedom) 1) that the stellar angular momentum 

points in the direction {6, 0) for the star IRS16 CC. The quantity p{x^, 1) = 1 — Erf (^x2/2) 
denotes the prob ability to observe by chance a value of or larger, even for a correct model 
JVao et al.ll2006[ ). 



For comparison, figure E] shows the distribution of reconstructed angular momentum 
directions for the same star. The maximum probability contour in figure M is identical in 
shape to the curve of the reconstructed angular momenta in figure [51 The 2;-prior determines 
the statistical weight of the curve in figure [51 

We define the deviation angle ip to be the minimum angle between the direction of the 
disk angular momentum and th e = contour. For srn all deviation angles, our definition 
is equivalent to the definition of iBeloborodov et al.l ( 20061). However, in our d efinition angles 
up to 180° are possible, while in the definition of IBeloborodov et al.l (l2006l ) the deviation 
angle is between zero and 90°. Figure [91 shows the distribution of the angular deviation ip 
of the disk stars from the common disk angular momentum direction for four MC generated 
disk of 1000 stars each, with 0°, 5°, 10° and 15° thickness. 
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Fig. 10. — (Left) cylindrical equal area projection of the distribution of significance in the sky for 
all 90 high quality WR/ stars with projected distances between 0.8" and 12". The disk positions 
of 
of 



Paumard et al 



Lu et al. 



1 200a ) are marked with full black circles and the position of the clockwise disk 
1 200A ) is marked by a broken black circle. There is a maximum s i gnific ance of 12.2a 
at {(j),0) = (262°, 51°), compatible with the clockwise system of \Paumard et al\ \200a ). The excess 
appears to have a narrow core and an extended tail to lower values of (p. Moreover, there is 
an extended excess of counter-clockwise orbits. (Right) Distribution of the negative cosine of the 
angular distance between the individual reconstructed angular momentum directions for each of the 
1000 generated z-values for all 90 WR/0 stars with respect to the observed excess center ((/>, ^) = 
(262°, 51°) (full blue line), the MC simulation expectation for 90 isotropic stars (red dotted line) 
and for a MC simulated disk of 35 stars with a two-dimensional sigma of 10°. 



5. Results 

5.1. The Clockwise System 

Figure [10] (left panel) shows a cylindrical equal area projection of the distribution of 
significances as a function of the sky coordinates cosOj vs. 0j for the 90 selected high 
quality WR/0 stars in the range < 12". There is a global maximum excess 

significance at {4>,6) = (262°, 51°) of 12.2cr. This corresponds to a reconstructed angular 
momentum density of 8.3 stars per 15° radius aperture. In the case of 90 isotropic stars 
the MC simulations predict an average value of 1.6 stars per 15° radius aperture, see figure 
E] (left panel) and an RMS of 0.55 stars per 15° radius aperture. The excess appears to 
have a narrow core. Figure [10] (right panel) shows the distribution of the negative cosine of 
the angular distance 6 between the individual reconstructed angular momentum directions 
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for each of the 1000 generated z- values for all 90 WR/0 stars with respect to the observed 
excess center. In addition it shows the respective distributions for MC simulated isotropic 
stars and a MC generated disk of stars with a two-dimensional sigma thickness of 10°. A flat 
distribution on a sphere is mapped to a flat distribution of the negative cosine of the angular 
difference. We determine the half width at half maximum (HWHM) of the excess as the 
angular difference S = 18°, for which the distribution of the negative cosine of the angular 
distance with respect to the excess center falls to half the maximum value. The width of the 
excess core is compatible with the width of a disk with a two-dimensional Gaussian sigma 
thickness (for a definition see equation [1]) of 10°. 

The peak angular momentum density for a MC simulated disk of 90 stars with a two- 
dimensional sigma thickness of 10° is 17.2 stars per 15° radius aperture (see figure [H] right 
panel). Comparing this value to the observed excess, we estimate that 35 stars contribute to 
this excess peak. Dividing the HWHM width of the excess peak by the square root of this 
number of stars, we estimate the statistical uncertainty of the peak position to be 3°. We 
re-ran the same analysis steps with Rq = 7.5 kpc, MsgrA* = 3.5 x lO^M© and Rq = 8.5 kpc, 
MsgrA* = 4.5xlO^M0. This yielded maximum offsets of 2° in and 6 and 0.5a in significance. 
Thus we estimate the systematic errors due to the uncertainties in Rq and Msgr a* to 2° in 
(f) and 6 and 0.5a in significance. Wi thin errors the excess peak position is compatible with 



the direction of the clockwise disk of iPaumard et al.l (120061 ). 



The excess has - in addition to the narrow core - an extended tail to lower values of 
Moreover, in figure [10] shows an extended U-shaped excess of counter-clock wise orbits 



including the direction of the counter-clockwise system of iPaumard et al.l (120061 ). 



5.1.1. Radial Dependence of the Excess Position 

In order to study a possible change in the excess features as a function of distance, 
we subdivided the sample of 90 WR/0 stars into three radial intervals: 32 stars in the 
interval 0.8"-3.5", 30 stars in the interval 3.5"-7" and 28 stars in the interval 7"-12". We 
chose the interval sizes such that each of the intervals contains approximately the same 
number of stars. Figure [11] shows cylindrical equal area projections of the significance sky 
distributions for these three intervals in projected distance to Sgr A*. In the inner interval 
there is a maximum excess significance of 13.9a at (0,^) = (256°, 54°), corresponding to a 
reconstructed angular momentum density of 5.5 stars per 15° radius aperture. In case of 
32 isotropic stars an average density of 0.63 stars and an RMS of 0.35 stars per 15° radius 
aperture are expected from MC simulations. The map shows a well-defined peak with a 
HWHM of 16°. It is compatible with the excess density of angular momenta from a MC 
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Fig. 11. — cylindrical equal area projections of the distributions of significance in the sky for three 
radial bins: (upper left panel) 32 WR/0 stars with projected distances in the bin 0.8" -3. 5" (upper 
right panel) 30 WR/0 stars in the bin 3.5" -7" and (lower panel) 28 WR/0 stars in the bin 7"- 
12". In the inner bin there is a max imum excess significa nce of 13. 9o- at {(l),9) = (256°, 54°), 
compatible with the clockwise system o uPaumard et all \200a ). The significance map in the middle 
interval shows two extended excesses, one for clockwise and one for counter-clockwise orbits. The 
clockwise excess has a local maximum significance ofhAa at {4), 9) = (262°, 48°) (compatible with the 
orientation of the clockwise system o uPaumard et al\ 11200a )) but a global maximum significance of 
5.9(7 at a clearly offset position: {(j),0) = (215°, 28°). The significance map in the outer bin shows a 
maximum excess significance of 11.5a at yet another position {(j),0) = (179°, 62°). The morphology 
of the excesses in the clockwise syste m may indicate a sm ooth transition of the excess center with 
projected radius. The disk positions iPaumard et al\\200a) are marked with black circles. 
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simulated disk of 25 stars with a 10° thickness. We thus estimate the error of the peak 
po sition to be 3.2 ° and fi nd the peak position to be compatible with the clockwise system 



of iPaumard et al.l (120061 ) . Only 4 out of the 32 stars in the inner radial interval are on 



counter-clockwise orbits. 

The significance map in the middle radial interval shows two extended excess structures, 
one for clockwise and one for counter-clockwise orbits. The clockwise excess structure has a 
local maximum significan ce of 5.4cr at ((t> , 6) = (262°, 48°) (compatible with the orientation 



of the clockwise system of iPaumard et al.l (120061 )) but a global maximum significance of 5.9a 



at a clearly offset position: ((^, ^) = (215°, 28°). 15 out of the 30 stars in the middle radial 
interval are on counter-clockwise orbits. An extended U-shaped excess structure at the 3- 
4(7 significance level is visible. For a more detailed discussion about this counter-clockwise 
excess see section E3i 

The significance map in the outer interval shows a maximum excess significance of 11.5cr 
at yet another position (0, 6) = (179°, 62°), corresponding to an angular momentum density 
of 4.6 stars per 15° radius aperture. For 28 isotropic stars, a density of 0.5 stars per 15° 
radius aperture is expected from MC simulations. The map shows a well-defined peak with 
a HWHM of 16°. It is compatible with the excess density of angular momenta from a MC 
simulated disk of 16 stars with a 10° two-dimensional sigma thickness. We estimate the 
error of the peak position to be 4°. 10 of the 28 stars in the outer radial interval are on 
counter-clockwise orbits. 

The angular distance between the significance peaks in the inner and outer intervals is 
(64 ±6)°. 

The significance map of the inner interval shows a significance below la at the maximum 
position of the outer interval {4>,6) = (179°, 62°). Moreover, the significance map of the 
outer interval shows a significance of only about 1.5o" at the maximum position of the inner 
interval {(f), 6) = (256°, 54°). We estimate the significance for a change of the maximum 
excess position as > lOcr. The morphology of the excesses in the clockwise system thus 
indicates that the excess center varies with projected radius. Table E] summarizes the excess 
positions and significances for the full sample of stars and the three radial intervals. 

We conclude that there is a significant change in the orientation of the clockwise system 
with projected distance. This change could refiect a large scale warp of a single disk, or 
the superposition of at least two disks or planar streamers located at different radii. The 
inner and outer radial interval can be well described by stars in 10° thick disks with a 
relative inclination of (64 ±6)°. The middle radial interval appears to represent a transition 
region. It also contains a significant fraction of stars on counter-clockwise orbits. In order to 
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Table 3: Parameters of the clockwise system. The given position error is the statistical error 
only. We estimate the systematic error as 2° in and 9. 



radial interval 


#stars 


max. significance 


max. 


position 


position error 


HWHM 


0.8" -12" 


90 


12.2a 




= (262°, 51°) 
= (98°, 129°) 


3° 


18° 


0.8" -3.5" 


32 


13.9(7 


(0,^) 

(0,0 


= (256°, 54°) 
= (104°, 126°) 


3.2° 


16° 


3. 5" -7" 


30 


5.9a 


(0,^) 

(0,0 


= (215°, 28°) 
= (145°, 152°) 


extended 


extended 


7" -12" 


28 


11.5a 


(0,^) 

(0,0 


= (181°, 62°) 
= (118°, 126°) 


3.8° 


16° 



explore this transition of the angular momentum direction of the clockwise system further, 
we ordered the 61 clockwise moving WR/0 stars by their projected distances to Sgr A*. 
We slid a window of width 19 stars over the ordered clockwise moving stars, resulting in 
42 groups of clockwise moving stars. For each group of stars we calculated the bin-by-bin 
sum of the x^(^, 0) sky plots for all the individual stars. We determined the sky position 
(^min, 0min) of the minimum of the summed sky plot. In the case of a warped system 
of stars, we would expect a smooth dependence of (6*111111, 0min) with distance. Isotropic 
stars would result in random fluctuations of (6*111111, ^min). In order not to be influenced by 
outliers, we iterate the determination of the minimum position. We calculated for all 19 
stars in the window the angular distance to the sky position (6*111111, 0min)- We excluded stars 
with a x^(6'min, 0min) > 100 and again computed the sky position with the minimum x^, 
(6'min,iter, 0min,iter)- This positiou is the average angular momentum direction of the group of 
stars. We note that the window selection on the projected distance may introduce biases in 
the average angular momentum position. Better selection criteria might be the 3D distance, 
the semi-major axis, or the total energy of the star in the potential well of the supermassive 
black hole. However, the use of these selection criteria introduces a model dependence as 
the 2;-position of the stars is unknown, not all stars are members of the clockwise system, 
and the candidate members have non-zero eccentricities. This will be the subject of future 
work. 

Figure [T2] (left) shows a cylindrical equal area projection of the local average angular 
momentum direction for the clockwise moving stars as a function of the average projected 
distance to Sgr A*. Figure [12] (right) shows the same data as gray scale points in orthographic 
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Fig. 12. — (Left) Cylindrical equal area projection of the local average stellar angular momentum 
direction for the clockwise stars as a function of the average projected distance. The points are 
correlated, see text. The average angular momeri t a for the innermost stars agree well with the 
orientation of the clock wise disk oflPaumard et al\ (200a ). shown by the black circle. The asterisk 
shows the Galactic pole i Reid & Brunnthale\ 200 A), th e diamond indicates the normal vector to the 



northern arm of the min isviral 'Jpaumard et all \2004) . the triangle indicates the normal vector to 



the bar of the minisyiral lLiszm2003i) . the square shows the rotation axis of the circum-nuclear disk 



(CND) backson et d\\l99,i ). (Right) Orthographic projection of the sky significance distribution 
of all WR/0 stars (see fiaure[W\). seen from (pQ = 226°, = 54° to highlight the clockwise system. 
The gray points show the average stellar angular momentum directions for the clockwise stars as 
a function of the average projected distance. The full black line shows a great circle between the 
angular momentum directions of the inner and outer borders of the clockwise system. The broken 
black line shows a fit with quadratic polynomials to 9j = 9j(R) and = (j)j{R). It describes the 
observed change of the angular momentum direction with radius better than the great-circle. 



projections centered on the clockwise disk together with the significance sky distribution from 
figure [TOl The average angular momentum direction of the clockwise stars is a function of 
the projected distance to Sgr A*. For the innermost stars, the average angu lar momentum 
direction is compatible with the clockwise system of lPaumard et al.l (l2006l ) but for higher 
projected distances the average angular momentum is clearly offset. This confirms the shift 
of the excess in the significance sky plots in figure [11] as a function of projected distance to 
the Galactic Center. The full black line shows a great circle between the angular momentum 
directions of the inner and outer borders of the clockwise system. The broken black line shows 
a fit with quadratic polynomials to Oj = 9j{R) and = (pj^R). It describes the observed 
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Fig. 13. — Local average stellar angular momentum direction for the clockwise stars as a function 
of the average projected distance: (Left) = 4>j{R) (Right) 9j = 6j{R). The J^2 points are 
correlated. They correspond to the ^2 positions of a window of width 19 stars which is slid over the 
61 clockwise moving WR/0 stars ordered by their projected distances to Sgr A*. 



change of the angular momentum direction with radius better than the great-circle. Figure 
Uni shows the local average stellar angular momentum direction for the clockwise stars as a 
function of the average projected distance: 0j = 4>j{R) (left panel) and 6j = Oj{R) (right 
panel). Within errors the angular momentum direction is compatible with being a smooth 
function of the average projected distance. We attribute the observed scatter to Poisson 
noise plus an additional local disk thickness. 

The most likely explanation for the observed change in the angular momentum direction 
of the clockwise system with distance to Sgr A* is a tilted warp in the clockwise disk (see the 
discussion section [6]), whose innermost part is the clockwise disk described in previous works. 
The angular difference between the innermost and outermost radii sampled is (64 ±6)°. The 
angular momentum direction as a function of the average projected distance does not follow 
a great circle, as expected for a simple warp. Instead, there is a significant offset from the 
great-circle, a tilt. 



Fig. 14. — Orthographic projection (seen from (po = 254°, = 54° J of the significance sky map 
in the interval of projected distances 0.8" -3. 5". It is overlaid with the 2 sigma contours (black 
lines) for the direction of the orbital angular momentum vect ors of the six early-typ e stars (S66, 
S67, S83, S87, S96 and S97) with 0.8" <R< 1.4" for which laHlessen et all kood ) were able to 
derive individual orbital solutions. All of these stars seem compatible with being members of the 
clockwise system. Still the orbital angular momenta of all of these stars are offset from the local 
angular momentum direction of the clockwise system at a confidence level beyond 90%. Th e whit e 
ellipse shows the 2 sigma contour of the c lockwise system as determined by \Paumard et al\ l200a ) 
and the brown one the 2 sigma contour of lLu et al\ ^2003 ). 



5.1.2. Inclinations to the Clockwise System 



Figure [14] shows an orthographic projection (seen from 0o = 256°, 9q = 54°) of the 
significance sky map in the interval of projected distances 0.8"-3.5" (upper left panel, figure 
[TTl) . It is overlaid with the 2a contours for the direction of the orbital angular momentum 
vector of the six early - type s tars (S66, S67, S83, S87, S96 and S97) with 0.8" <R< 1.4" for 
which iGillessen et al.l (l2009l ) were able to derive orbital solutions. All of these stars appear 
to be compatible with being members of the clockwise system. Still the orbital angular 
momenta of all of these stars are offset from the local angular momentum direction of the 
clockwise system at a confidence level be yond 90%. The white ellipse shows the 2a contour 
of the clockwise system as determined by lPaumard et al.l (120061 ). The mean angular distance 
between the angular momenta of the six stars and the disk angular momentum direction is 
< %jj >= 14.5° and the mean of the angular distances squared is < ip'^ >i/2= 15.4°. The 
distribution of angular offsets can be fit with a two-dimensional Gaussian (see equation [1]) 
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Fig. 15. — Distribution of the reconstructed angular difference ijj for all the 90 WR/0 stars from 
the local average angular momentum direction of the clockwise system (blue histogram). The peak 
is well described by the expected distribution for a disk with a 10° two-dimensional Gaussian sigma 
thickness and 46 disk members. In addition, there is a long tail to large angular distances towards 
to clockwise system. 

with = (11 ± 2)°. 

In a next step we want to determine the angular offsets of all WR/0 stars from the 
local angular momentum direction of the clockwise system. We interpolated the points 
shown in figure [12] (left) with a third degree spline function in projected distance to obtain 
the direction of the clockwise system as a function of average projected distance. Figure 
[TSl shows the distribution of the reconstructed angular difference ip for all the 90 WR/0 
stars from the local average angular momentum direction of the clockwise system (blue 
histogram). The distribution has a peak at small inclinations, which is well described by 
the expected distribution for a disk with a 10° two-dimensional Gaussian sigma thickness 
and 46 disk members. In addition, there is a long tail to large angular offsets from the 
clockwise system. 11 out of the 61 clockwise stars have inclinations larger than 30° from the 
angular momentum direction of the clockwise system. The clustering between 80° and 100° 
inclinations is due to the counter-clockwise rotating stars. 

Within errors the inclinations of the stars to the clockwise system are independent of 
distance to Sgr A*. 
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Fig. 16. — Surface number density of the 30 candidate stars in the clockwise system (blue filled 
circles) as a function of the three-dimensional distance from Sgr A* , which have a minimum 
angular distance below 10° from the (local) average an gular momentum direc tion of the clockwise 
system. The dashed line shows the r~^'^ power-law of lPaumard et ali /(200a ). The full blue line 
shows the best fit power-law to the cock-wise disk: T,{rdisk) oc r~ 



-1.95±0.25 



5.1.3. Surface Density 

Both the calculation of the surface density of the clockwise system as well as the calcu- 
lation of the stellar orbital elements requires the determination of the unknown ^-coordinate. 
In order to obtain the least biased results, we apply a hard criterion for the choice of can- 
didate stars for the clockwise system. We only consider the 30 of the 61 clockwise moving 
WR/0 stars which have angular distances below 10° from the local angular momentum di- 
rection of the clockwise system. The fraction of stars which have angular distances below 10° 
from the local angular momentum direction of the clockwise system of all clockwise moving 
stars is independent of distance to Sgr A* within errors. For each star we generate 1000 z 
positions taking into account the probability for the stellar angular momentum direction 
(see e.g. figure [8]). Moreover we apply the prior that the local inclinations to the clockwise 
system follow a two-dimensional Gaussian distribution with a = 10°. 

We compute the surface number density of observed candidate stars of the clockwise 
system as: S(rdisk) = AiVstars(?^disk)/AAdisk(?^disk) with AAdisk(?^disk) = 27rrdiskArdisk, where 
Tdisk is the distance to Sgr A* measured in the local angular momentum direction of the 
clockwise system. Figure [12] shows the surface number density of candidate stars in the 
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clockwise system as a function of the three-dimensional distance from Sgr A* in the local 
angular momentum direction of the clockwise system. The best fit power-law to the clockwise 
system is S(rdisk) oc '"disif estimate the systematic error of the power-law slope due 

to uncertainties in Rq an d Msgr a* to 0.2. This res ult is con i patibl e with the value found 
bv IPaumard ^Tdl J2006h . S(rdisk) oc r-f^^o^ and ILu et all J2008h . S(rdisk) oc r-^^^^-^^- 
Both groups did not take into account the change of the angular momentum direction with 
distance from Sgr A*. 



5.1.4- Reconstruction of Orbital Elements 

The analysis above was, to lowest order, independent of the enclosed mass. The potential 
only infiuenced the range of the generated z-values due to the assumption of bound orbits, 
see section 14.1.11 On the contrary, the determination of stellar orbital elements depends on 
the enclosed mass of the orbit, which is given by the sum of the mass of Sgr A* plus the 
mass of the cluster of late- type stars. We used the following parameter ization of the enclo sed 



stellar mass M(rstar) as a function of the 3D distance rstar to Sgr A* (ITrippe et al.ll2008l ): 



fi / 2.1 X lO^Mopc-3 , 



(r/8.9")2 

The presence of an ex tended mass indu ces retrograde pericenter shifts which result in open 



rosetta shaped orbits (lAlexanderl l2005l ) . We approximate such an orbit locally by a Kepler 
ellipse and use the enclosed mass in a sphere given by the instantaneous position of the star 
to compute the orbital elements. 

Figure [T7| shows the distributions of the orbital elements a, e, u, t for the 30 candidate 
stars of the clockwise system. The distribution of semi major axis a (upper left panel of figure 
[T71) is compatible with an power law. The distribution of the reconstructed eccentricities 
(upper right panel of figure [T7j) has a mean of 0.51 and an RMS of 0.27. The eccentricity 
distribution shows a two-peak structure: a broad distribution between and 0.8 and five 
stars between 0.9 and 1. Figure [TS] shows the eccentricity of the 30 candidate stars as a 
function of projected distance to Sgr A*. The five stars with high eccentricities have very 
small RMS values for their reconstructed eccentricities. They could only be disk members in 
a very small region of their parameter space, limited to nearly radial orbits. This fine-tuning 
may be a sign that these five stars do not dynamically belong to the clockwise system, but 
rather are by chance compatible with being candidates. Still, we cannot exclude them at 
this point. Without these five stars, the mean eccentricity is 0.42 ± 0.05 and the RMS is 
0.20 ± 0.03. We note that the eccentricity is a purely positive quantity. For bound orbits it 
is less than one. Therefore the eccentricities do not follow a Gaussian distribution. We use 
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Fig. 17. — Distributions of the orbital elements for the 30 clockwise moving WR/0 stars in the 
radial bin 0.8"-12", which have a minimum angular distance below 10° from the (local) average 
angular momentum direction of the clockwise system. (Upper left) dNstars/da, the full line shows an 
power law, (upper right) dA^stars/de, (lower lef t) dA^stars/do;, (lower right) d iYstors/dgfAffrstar), 
g = 2t for Vz > and g = 2(1 — r) for Vz < 0, see Weloborodov & Levin 1(2004) ■ 



the observed mean and RMS eccentricities to quantify the first two moments of the a priory 
unknown eccentricity distribution. Projection effects and the prior of a clockwise system with 
a Gaussian distribution of inchnations bias the distributions of reconstructed eccentricities to 
larger values. We have applied our reconstruction method to MC simulated stars in a disk 
with a 10° two-dimensional Gaussian sigma thickness with a flat eccentricity distribution 
(mean eccentricity of 0.5 and RMS of 0.28). The reconstructed eccentricity distribution has 
a mean of 0.55 and an RMS of 0.26. We estimate the bias of the average reconstructed 
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Fig. 18. — Reconstructed eccentricity as a function of projected distance for the 30 clockwise 
moving WR/0 stars (blue points), which have a minimum angular distance below 10° from the 
(local) average angular momentum direction of the clockwise system. Red ci rcles show the six early - 
type stars (S66, S67, S83, S87, S96 and S97) with 0.8" <R< 1.4" for whic Aoillessen et al\ fmM) 
were able to derive individual orbital solutions. Error bars denote the RMS of the reconstructed 
eccentricities. 



eccentricity to 0.05 and in addition quote a systematic error of 0.05. We estimated the 
systematic errors due to uncertainties in Rq and Msgr a* to be as small as 0.02. The mean 
eccentricity of the candidate stars of the clockwise system is 0.37 ± 0.05stat. ± 0.05s 



syst. 



Moreover, figure [T8] includ es the six early- t ype st ars (S66, S67, S83, S87, S96 and S97) 
with 0.8" < R< 1.4" for which iGillessen et al.l ( 120091 ) were able to derive individual orbital 
solutions (see figure [H]). Of these six so-called S-stars only two (S87 and S96) fulfill our 
strict selection criteria for disk candidate stars. The reconstructed eccentricities of S87 and 



596 in this work agree within errors with the values of IGillessen et al.l (120091 ). S67, S83, and 

597 have inclinations of more than 10° to the disk angular momentum direction and S66 is 
a B dwarf {mK = 14.8). The six S-stars have a mean eccentricity of 0.36 and an RMS of 
0.20, giving an error of the mean of 0.09. In contrast the so-called S-stars have a different 
eccentricity distribution: dA^'stars/de oc e^-eio.g (IGillessen et al.ll2009l ). We con abine our result 



of the mean eccentricity with the mean eccentricity of the six stars from IGillessen et al. 
(I2OO9I ) to a weighted average of 0.36 ± 0.06. 



Figure [TTT lower left panel) shows the distribution of the reconstructed arguments of 
periapsis of the orbits u. For uniformly populated disks and a uniform azimuthal exposure 
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the distribution of uj is expected to be flat. The observed distribution is compatible with a 
flat distribution within errors. 

Figure [TTIflower right panel) shows the distribution of the time separating the stars from 
their nearest passage of their pericenter normalized to their half-periods 5'(M(rstar))- For the 
true enclosed mass and a random snapshot time, the expected g obeys Poisson statistics: it 
has a flat probability distribution between and 1 with the mean expectation value ^ = 0.5 
and the standard deviation A q = 12~^/^ 0.29. q — » for small assumed masses an d 'g ^ 1 



for too large assumed masses (iBeloborodov fc Levin! I2004J : iBeloborodov et al.l l2006l ) . In our 
data we determine ^(M(rstar)) = 0.52 ± 0.05 and /^g{M{r^t^r)) = 0.27 ± 0.04. In case we 
adopt Rq = 7.5 kpc and Msgr a* = 3.5 x IO^Mq we get g = 0.48±0.05 and for the parameters 
Ro = 8.5 kpc and Mggr a* = 4.5 x lO^Mg we get g = 0.53 ± 0.05. 



5.2. Hertzsprung-Russell Diagram 



Figure [19] shows the location of the 0-type super giants, giants and dwarfs and B super 
giants in a Hertzsprung-Russell (HR) diagram. Such a diagram can be used to estimate the 
age of the population by comparing of the stellar positions to theoretical isochrones. Such 
isochrones are displayed as solid grey lines. Our HR diagram excludes Wolf-Rayet stars 
since isochrones are ill-defined at the position occupied by these stars. For each star we 
assigned effective temperatur es derived from spectral types us i ng th e effective temperature 
scales of iMartins et al.l (120051 ) for O stars and lCrowther et al.l (120061 ) for B supergiants. We 
es timated the luminositie s from absolute K-band magnitudes and the bolometric corrections 
of iMartins fc Plea (120061). We assumed a distance of (8 .0 ± 0.5) kpc and an extinction of 
Ak = 2.8 ± 0.5 (IScoville et al.ll2003l : ISchodel et al.l 120071 ). The errors in Tgfr correspond to 
uncertainties in spectral types. These er r ors pr opagate into the bolometric corrections, see 
equation 3 and table 1 of iMartins &: Plea (120061 ) . and lead to the uncertainty in luminosities. 
In figure [191 we differentiated four subgroups of stars: the candidate stars of the clockwise 
system with a maximum angular offset of 10° from its local angular momentum direction, 
the clockwise rotating stars which have offsets between 10° and 30° and the one with offsets 
larger than 30° from the clockwise system, as well as the counter-clockwise stars. Stars with 
log L/L0 < 5.0 and log Teg < 4.45 have a significant probability to be younger than 8 Myr. 
Indeed, these stars have the same probability to have logTefj = 30000 K (corresponding to 
an age of 4-8 Myr) or logTgfj = 20000 K (corresponding to older ages). This is due to the 
absence of Tes indicator in the K-band spectra of stars with 20000K < Teg < 30000 K. 
Besides, the small number of red supergia nts argues against an age much larger than 10 Myr 
JPaumard et al.lbood : Bavies et al.ll2008h . 
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Fig. 19. — Location of the 0-type super giants, giants and dwarfs and B super giants in a 
Hertzsprung-Russel diagram. Red triangles denote the candidate stars with a maximum angular 
offset of 10° from the local angular momentum direction of the clockwise system, the orange pen- 
tagons show the clockwise stars with off sets between 10° and 30° and the blue points show the 
clockwise rotating stars which have a minimum offset of 30° from the clockwise system. The green 
squares represent the counter-clockwise stars. The numbers indicate the initial mass of the stars 
on the tracks starting from the zero age main sequence (solid black line). The two thick gray lines 
show two isochrones f or 4 and 8 Myr. The trac ks (various dotted/dashed lines) are based on Geneva 
models with rotation iMeunet & Maeden\200a ). 



The main conclusion is that within errors these four groups o f stars are coeyal wit h 
an age between 4 and 8 Myrs, comparable to the earlier studies by iPaumard et al.l (120061 ). 
There is no variation of average stellar age with projected distance to Sgr A*. 
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Fig. 20. — (Left) cylindrical equal area projection of the density of reconstructed angular momenta 
in a fixed aperture of 25° radius for 30 WR/0 stars in the interval of projected distances to the 
Galactic Center 3. 5" -7". There is an extended U-shaped excess of counter-clockwise orbits with 
a global maximum angular momentum density of 3.8 stars per 25° radius aper t ure a t {<j),6) = 
(200°, 142 °), near the position of the counter-clockwise system of \Paumard et al\ l200a ). The disk 
positions iPaumard et al\\200a ) are marked with black circles. (Right) Probability that an isotropic 
distribution of stars yields in at least one aperture of 25° radius a certain density of reconstructed 
angular momenta. A maximum angular momentum density of 3.8 stars per 25° radius aperture 
corresponds to a probability of 2%. 



5.3. The Counter-Clockwise System 

The exi s tence of a co unter-clockwise disk of early-type stars has been called into question 



bv iLu et al.l (120061 . 120081 ). It is undisputed that about 30% of aU WR/0 stars (29/90) are 
on counter-clockwise orbits. So the question is not whether such stars exist but whether 
their kinematic distribution is indicative of a second, well defined disk structure. The upper 
right panel (interval of projected distances 3.5"-7") of figure [TT] shows an extended U-shaped 
excess in the distribution of counter-clockwise stellar angular momenta. We calculated the 
significance for a sky density of angular momenta averaged over an aperture of 15°, chosen 
to maximize the significance of a moderately extended disk of 10° two-dimensional sigma 
extension. Figure [20] (left panel) shows for the same stars as in the upper right panel 
(interval of projected distances 3.5"-7") of figure [TT] the density of reconstructed angular 
momenta in a fixed aperture of 25° radius, better suited for more extended features. There 
is an extended U-shaped excess of counter-clockwise orbits with a global maximum angular 
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Fig. 21. — Cylindrical equal area projections of the density of reconstructed angular momenta in 
an aperture of 25° radius (left) average for a MC simulated disk of many stars, sc aled to 15 stars 



The d isk angular momentum points in the direction of the counter-clockwise disk o uPaumard et al 



1 200a) . (right) 15 MC simulated stars with the same azimuthal distribution as the observed stars. 



All 15 stars are members of a disk in the direction of the counter-clockwise disk of lPaumard et al 
1 (200(1 ). 



momentum density of 3.8 stars per 25° ra dius aperture at (4), 9) = (200°, 142°), near the 
position of the counter-clockwise system of lPaumard et ahl (120061 ). In case of 15 isotropic 
distributed counter-clockwise stars we expect a mean density of 1.5 stars and an RMS of 0.52 
stars per 25° radius aperture centered at (0, 9) = (200°, 142°), corresponding to an excess of 
4.5 times the RMS. In figure [2D] (left panel) we did not use an y additional z-position prio r 
for the IRS 13E stars. If IRS 13E is a star cluster at ;z ~ 7" JPaumard et al.lbooi l2006h . 
the stars IRS 13E2 and E4 would not be candidate members of a counter-clockwise system. 
In this case, the maximum angular momentum density of counter-clockwise orbits would 
slightly be reduced to 3.7 stars per 25° radius aperture at the same position as given above. 

There are 15 stars with counter-clockwise orbits located at projected distances between 
3.5" and 7". The location of the excess is not known a priori. What is the probability to 
find such an excess in an isotropic distribution of stars when searching all possible aperture 
directions? We computed angular momentum density maps (25° radius aperture) for many 
groups of 15 simulated isotropic stars on counter-clockwise orbits. Figure [20] (right panel) 
shows the probability that an isotropic distribution of stars yields, in at least one aperture 
of 25° radius, a certain density of reconstructed angular momenta. A maximum angular 
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Fig. 22. — (Left) Cylindrical equal area projection of the local average stellar angular momentum 
direction for the counter-clockwise stars as a function of the average projected distance. The points 
are correlated, see text. The average angular momenta for the innermost 16 stars (corresponding 
to the 7 points with sm allest average p r ojecte d distance) agree well with the orientations of the 
counter-clock wise disk of\Paumard et al\ Il200a ). shown by the black circle. The asterisk shows the 



Galactic pole iReid & Brunnthalei 



200 A), the diamond indicates the normal vector to the northern 



arm of the min ispiral IPaumard et al. 2004 )- the triangle indicates the norm al vector to the bar o f 



the minispiral l(Liszi\2003 ). the square shows the rotation axis of the CND ^Jackson et al\\l993i ). 



(Right) Orthographic projection of the density of reconstructed angular momenta in an aperture of 
25° radius for the 30 WR/0 stars in the interval of projected distances 3.5"-7" (see figure 
seen from (pQ = 160°, Oq = 142° to highlight the counter-clockwise stars. The gray points show 
the average stellar angular momentum direction for all counter-clockwise stars as a function of the 
average projected distance. 



momentum density of 3.8 stars per 25° radius aperture corresponds to a probability of 2%. 
Taking the whole observed U-shaped excess structure into account the probability is even 
lower. Hence we conservatively exclude at the 98% confidence level that the observed excess 
of counter-clockwise stars is due to a fluctuation of isotropic stars. 

In addition to the global maximum angular momentum density, the significance sky 
distribution (left panel of figure [201) shows an extended U-like shape with local maxima up 
to 3.7 stars per 25° radius aperture. This U-like shape may be due to 

• projection effects in the reconstruction of nearly face-on orbits (cos^ ~ —1) 
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Fig. 23. — Local average stellar angular momentum direction for the counter-clockwise stars as a 
function of the average projected distance: (Left) = (pj{R) (Right) 9j = 9j{R). The 20 points 
are correlated. They correspond to the 20 positions of a window of width 10 stars which is slid over 
the 29 clockwise moving WR/0 stars ordered by their projected distances to Sgr A* 



a small number of stars with a possible non-isotropic distribution on the sky 

an intrinsic property of the counter-clockwise structure of WR/0 stars like a large 
thickness, a warp, a state of dissolution or the presence of separate streamers etc. 



Figure [^ITleft) shows the average density of reconstructed angular momenta in an aper- 
ture of 25° radius for a MC simulated disk (10° thickness) of many stars, scaled to 15 
st ars. The disk angular momentum points in the direction of the counter-clockwise disk 
of iPaumard et al.l (120061 ). Although the simulated disk is nearly face-on (cos 6' ~ —1) the 
angular momentum direction of the disk is not biased by the reconstruction procedure. Fig- 
ure [2T](right) shows one representation of 15 MC simulated stars with the same azimuthal 
distribution as the observed counter-clockwise stars. All 15 stars are members o f a disk (10° 
thickness) in the direction of the counter-clockwise disk of iPaumard et al.l (120061 ). The excess 
has the same U-like structure as observed in the data with similar values for the angular 
momentum density. Hence the combination of a small number of observed stars, which are 
members of a 10° thick disk, at their observed positions on the sky can — but does not 
necessarily has to — produce the U-like shape of reconstructed angular momenta. 



In an analogous study to the clockwise system, we investigated the average angular 
momentum of the counter-clockwise stars as a function of projected distance, keeping in 
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mind that the small number of stars limits the significance of the result. We ordered the 
29 counter-clockwise moving WR/0 stars by their projected distances to Sgr A*. We slid 
a window of 10 stars width over the ordered counter-clockwise moving stars and computed 
the average angular momentum direction for these stars. Figure [22] (left) shows a cylindrical 
equal area projection of the local average angular momentum direction for the counter- 
clockwise moving stars as a function of the average projected distance to Sgr A*. Figure 
|22] (right) shows the same data as gray scale points in an orthographic projection centered 
on the counter-clockwise excess together with the density of reconstructed angular momenta 
in an aperture of 25° radius. The average angular momentum direction of the innermost 
16 counter-clockwise stars (corresponding to the 7 windows of 10 stars each with smallest 
average project e d dist ances) clusters around the position of the counter-clockwise system of 



Paumard et al.l (120061 ). The average angular momentum direction of the counter-clockwise 
stars may be a function of the projected distance. Figure [23] shows the local average stellar 
angular momentum direction (0j = 0j(-R) and 6j = 6j{R)) for the counter-clockwise stars as 
a function of the average projected distance. Due to the small number of counter-clockwise 
stars, the spread is large and the significance of a possible shift of the average angular 
momentum direction is marginal. Still the plot 0j = (f)j{R) suggests an evolution of the 
(p angle of the angular momentum vector with distance to Sgr A*. 10 of the 29 counter- 
clockwise stars have inclinations larger than 30° with respect to the local angular momentum 
direction of the counter-clockwise system. 



Our orbita l analysis is in good agreement with the results of iGenzel et al.l (120031 ) and 



Paumard et al.l (120061 ) which were only based on velocity vectors. They found that there is 
probably a significant second, counter-clockwise system, at large angles with respect to the 
main clockwise system. The significance of this second system is moderate. The observed 
maximum angular momentum density occurs with a 2% chance in an isotropic distribution 
given a large number of tries (45 sky bin for counter-clockwise orbits cos^ < 0). We cannot 
exclude that the observed counter clockwise structure is due two separate streamers etc. 
However, the significance of the counter-clockwise system is about as large as can be rea- 
sonably expected on the basis of the small numbers of stars contained in it, as well as its 
observed azimuthal and broader angular distribution of stars. 



5.4. Summary and Comparison to Previous Results 



In this work we have increased the data set by 27 new reliably measured WR/0 stars in 
the central 12" of our Galaxy, to a total of 90, and reduced the prop er motion uncer t aintie s 
by a factor of four compared to the data set previously analyzed by [Paumard et al.l (|2006[ ). 
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The p revious analysis by iLevin fc Beloborodovl (120031 ) ; iGenzel et al.l (120031 ) ; iPaumard et al. 
( 120061 ) only included the 3D velocity information of the early-type stars to compute a for 
stellar disks. Here we have i n addi tion taken into account the 2D positions on the sky. The 
results from iPaumard et al.l (120061 ) were based on 59 quality 1 and 2 early-type stars with 
projected distances between 0.8" and 8". 

We have shown in this paper that the angular momentum direction of the clockwise 
system is a function of the projected distance from Sgr A*. Its innermost edge is statisti- 
cally compatibl e with the angular momentum direction of the clockwise disk determined by 



Paumard et al.l (120061 ). There does not appear to be any significant difference between an 
analysis with all velocitie s and spatial co o rdina tes and an analysis based on velocities only. 
In a reanalysis of the full IPaumard et al.l (120061 ) data set (all quality 2 stars with projected 
distances between 0.8" and 12") we se e already a hint of th e angular momentum direction 
change, albeit at lower significance. As IPaumard et al.l (120061 ) restricted their analysis to the 
range of projected distances between 0.8" and 8", they did not see any significant change of 
the average angular momentum. 

The surface number densit y of observed stars is within errors compatible with the one 
given by IPaumard et al.l (120061 ). Also the distribution of inclinations of the stellar orbits 
to the local angular momentum di rection of the clockwis e system agrees well with the disk 
thickness of (14 ±4)° determined by IPaumard et al.l (120061 ). The average star in the clockwise 
system is not on a circular orbit. Our analysis indicat es somewhat larger ec centricities of 
the candidate stars of the clockwise s ystem compared to IPaumard et al.l (120061 ). Our results 



are compatible with the analysis of iBeloborodov et al. 



early-type stars with 0.8" < R < 1.4" by iGillessen et al. 




and the results for the six 



Paumard et al.l (120061 ) used a uniform z-prior to compute the eccentricities and noted 



a bias to too larges value s, which they t r ied to correct by applying a weight of 1/(1 — 
exp(— e/0.35)). Moreover, IPaumard et al.l (120061 ) used a smaller mass of Sgr A* of 3.6 x 
10^ Mq and did not take the mass of the late-type cluster into account. A smaller mass 
favors smaller eccentricities. Firmer conclusions on the eccentricity distribution will require 
a more thorough analysis of the systematic effects in the eccentricity reconstruction. 

In addition, we confirm a non-isotropic distribution of the counter-clockwise stars at a 
confidence level beyond 98%. The modest significance of this counter-clockwise structure 
is a combination of the small number of stars (Poisson noise) and the fact that it appears 
to be less well defined. The counter-clockwise system may be a second disrupted disk or 
streamer of stars. For simplicity of analysis and visualization we have described the stellar 
structures in the central parsec in our discussion above in terms of fiat disks or rings. We 
call deviations from such simple flat structures 'disk warps' or 'dissolving/disrupted disks'. 
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This approach is based on our aim of distilhng the simplest physical model required and 
justified by the data. The reality is undoubtedly more complex. The stellar distribution in 
the Galactic Center may consist of a connected but complex system of separate streamers 
at different radii and different orientations. 



Lu et al.l (120081 ) independently performed a similar analysis searching for disk features 



amongst the early-type stars in the Galactic Center. They use two different data sets: a 
primary set of 32 stars at the projected distance bin 0.8"-3.5", for which they measured 
p ositions, proper motion velocities and acceleration upper limits but use the radial velocities 
of iPaumard et al.l (120061 ). In this radial distance bin they find that most stars are members 
of a fairly thin (RMS dispersion angle of (7 ± 2)°) clockwise disk at ^2 = (100 ± 3)° and 
i = (115 =b 3)°. These numbers are in reasonable agreement with ours (above), as well 



as 



Paumard et al 



(120061 ). The slight offset in inclination may be due to different priors 
used for the distribution of the unknown z-coordinate of the young stars. ILu et al.l (120081 ) 
do not use the volume density of the early-type stars to generate the 2;-coordinate (see 
section I4.1.ip but rather use a uniform acceleratio n prior. Moreov er, the dispersion angle 
depends on the definition of candidate disk stars. ILu et al.l (120081) show that some of the 
candidate star s of the clockwis e system have eccentricities larger than 0.2 in agreement with 
our analysis. ILu et al.l (120081 ) also analyze a second (extended) data set of 73 stars: the 
32 stars of thei r prim ary data set and in addition 41 quality 1 and 2 early-type stars from 
Paumard et al.l ( 120061 ) at projected distances between 3.5" and 12". They concluded that 
the number probability of candidate stars being members of the clockwise disk decreases 
with increasing projected distance from Sgr A*. Still they do not find this decrease to be 
significant at a level beyond 3a. Moreover, they do not find the pre-trial significance of the 
counter-clockwise system above la. There are five main differences, why our analy sis yields 



a high er significance of the counter-clockwise system compared to the analysis of iLu et al. 



the counter-clockwise disk has very few candi date merabers w ith projected distances 
below 3.5", the limit of the primary sample of iLu et al.l (|2008l ) 



analyzing the extended sample, ILu et al.l (120081 ) do not make cuts on the projected 
distance 

in this work we have reduced the proper motion u ncert ainties by a factor of four com- 



pared to the values published by iPaumard et al.l (120061 ) , which were used by ILu et al. 
( l2008h 



we added 27 new reliably measured WR/ O stars to the data set, 15 of them on counter- 
clockwise orbits 
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we required a firm identification as a WR/0 sta r with a maxim um radial velocity 



uncertainty of 100 km/s. In their extende d aamp/e iLu et a. 



1120081) included 3 B stars 



and one O star affected by crowding from lPaumard et al.l (120061 ) . which we do not use 
in our analysis. 



Given the larger number o f stars an d sma ller errors in our data set for projected distances 
beyond 3.5", the results of ILu et al.l (120081 ) are not in contradiction with our results. 



6. Discussion 



Our an alysis confirms the existence of a relatively thin main disk-like structur e of clock- 



wise stars ( iLevin &: Beloborodovl l2003l : iPaumard et al.l l2006l : iTanner et al.l l2006l : ILu et al. 



20081 ). This structure exhibits a significant change of its rotation axis at different radii. 
The main disk may be strongly warped, or consists of several, closely related streamers with 
radially varying orbital planes. We also confirm at th e 98% confiden ce level the existence 



of a coherent structure of the counter-clockwise stars (IPaumard et al.ll2006l ). This second 
system may be a set of streamers at large angles relative to the primary clockwise system, 
or perhaps a second, albeit more disrupted, disk. 

In the following we discuss the possible nature and origin of these disk-like structures in 
the framework of different models. Any such model can be tested against several observables: 



• number of disks, fractions of disk candidates and random stars, disk thickness 

• disk orientation, disk warps 

• radial distribution of the early-type stars 

• eccentricity distribution of disk candidate stars 

• mass function. 



We will shortly summarize our data on these five observables, and then discuss them in 
view of both the inf ailing cluster and the in situ star formation models, described below. 



Fig. 24. — Three-dimensional visualization of a warped and tilted clockwise disk. It is shown 
assuming locally flat circular orbits (a = 1"-10"). The angular momentum vector is a smooth 
function of semi-major axis, the broken line in figure{T^ (right panel). For comparison, the flat, 
non-warped ring in the range a= 3.5"~7" at large inclination to the clockwise warped and tilted disk 
indicates the average angular momentum direction of the innermost 16 counter-clockwise stars. 



6.1. Observables 

6.1.1. Number of Stellar Structures and their Thickness 

55% of the observed WR/0 stars in the central 0.5 pc of the Galactic Center are 
candidate members of a warped clockwise system and 20% of the stars are candidate members 
of a counter-clockwise system at an inclination of about 100° with respect to the clockwise 
system. Both systems are not locally fiat, the inclinations to the clockwise system are 
compatible with a two-dimensional Gaussian distribution with a sigma of 10°. Moreover, 
25% of the WR/0 stars have angular separations from the disks too large to be considered 
candidate members. However, there is no further noticeable difference between the disk star 
candidates and the apparently random stars given their similar ages, see section 15.21 The 
data thus suggest the existence of two main structures; a main clockwise relatively thin 
disk-like structure and a somewhat thicker less massive counter-clockwise structure. 
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6.1.2. Warps and Orientation 

The rotation axis of the observed clockwise system changes as a function of the projected 
distance from Sgr A*. Its direction changes by Ag = 26° and = 41° between the inner 
region and the intermediate one; and then changes again by Ag = —34° and = 36° in the 
transition to the outer region (see table [3D • The angular momentum direction as a function of 
the average projected distance does not follow a great circle as expected for a simple warp. 
Instead, there is a significant offset from the great-circle, a tilt. The radial dependence 
of 6 J and 0j can reasonably well be fitted by quadratic polynomials, see figure [12] (right 
panel). The rotation axis of the counter-clockwise structure keeps a constant inclination 
(6) of approximately 100° from the main clockwise disk but has a broad distribution in 
throughout the intermediate region where it is observed. 

Figure [23] shows a three-dimensional visualization of the warped and tilted clockwise 
disk. It is shown assuming locally flat circular orbits (a = 1"-10"). The angular momentum 
vector is a smooth function of semi-major axis, given by the quadratic fits to the direction 
of the angular momentum of the clock- wise system as a function of distance to Sgr A*. For 
comparison the flat, non-warped ring in the range a= 3.5"-7" at large inclination to the 
clockwise warped and tilted disk indicates the average angular momentum direction of the 
counter-clockwise stars with projected distances between 3.5" and 7". 



6.1.3. Eccentricity Distribution of Disk Candidate Stars 

The main clockwise structure is very coherent and is most likely a single cohesive disk of 
stars. The eccentricity of most likely disk candidate stars can therefore be estimated under 
the assumption that they indeed belong to such a disk. The data suggest non-zero, moderate 
eccentricities. The a priory unknown eccentricity distribution has a mean of 0.36 ± 0.06 and 
an RMS of 0.2. We can exclude at a confidence level beyond 5 sigma that all disk stars are 
on circular orbits. There are four candidate disk stars with eccentricities above 0.9, which 
may belong to a different population of stars (e.g. they are no disk stars). 



6.1.4- Mass Function (MF) and Radial Distribution 



Paumard et al.l (120061 ) showed that the MF of the disk sta rs is significantly flatter than a 



Salpeter MF. Also the observed X-ray emission constraints of iNayakshin fc SunyaevI (120051 ) 
suggest a top heavy MF. The surface number density of the disk stars is a steep function of 
distance from Sgr A*: S(rdisk) oc The current limited data on the mass function 
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of the WR/0 stars suggest the same distribution of the MF throughout the main clockwise 
system at different radii. It also seems to be similar to the MF of stars in the counter- 
clockwise system. However, to better constrain the MF of these stars, the identification of B 
dwarfs throughout the regions where the young stars are observed is required, in addition to 
the dynamical analysis of their kinematics. Moreover, a thorough completeness correction 
is required for a better understating of the MF. These caveats are beyond t he scope of this 
paper and will be addressed in a forthcoming paper by iMartins et al.l (120091 ). 



6.2. Models for the origin of the young GC WR/O stars 



The two main scenarios discussed in the literature for the origin of the disk-like structure 
of young stars in the GC are the inf ailing cluster model and the in situ star formation model. 
In the infalling cluster scenario, a cluster of young stars formed far from the GC (a few 
parsecs) and in-spirals to Sgr A*. During the in-spiral the stars are stripped from the young 
cluster, forming a disk-like structure. In a mass segregated cluster the most massive stars 
are positioned in the cluster core and are the last to be stripped at the innermost regions 
of the GC, close to the MBH. The final stripped stars are suggested to be the young O and 
WR stars observed in the GC. 

In the in-situ model, a clump or clumps of gas fall into the GC, where they form a 
disk-like structure. The gaseous disk then fragments to form new low mass stars in a stellar 
disk, that grow, through accretion, into the more massive stars currently observed in the GC 
flhevin fc Beloborodoyll2008l: lGenzel et ahlbooal: lGoodmanl[2008l: iMilosav^^ fc Loebl 2004': 

20081 : 



J 



Navakshin fc Cuadra 



Hobbs fc Nayakshin 



2005 



20081 



Paumard et al 



200 



Bonnell fc Rice 12008 



Mapelli et a] 



Any valid model can be constrained by the observations and should be able to explain 
the observed data described above. In the following we confront the current models with the 
data. 



6.2.1. In situ formation 



The tidal forces near the super-massive black hole Sgr A* are so strong that a cloud can- 
not be gravitationally bound unless its density exceeds 6-lO^MsgrA*/ (4-lO^M0) (i?/7" )~^cm~^, 
where R is th e distance to Sgr A* , orders of magnitude denser than currently observed 
( iMorrisI Il993l ). The tidal shear can be overcome if the mass accretion was large enough 
at some point in the past causing a gravitationally unstable massive (~ IO^Mq) disk to 
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form. The stars formed directly out of the fragmenting disk. Given the fact that the star 
formation event 6 Myrs ago seems to be i solated, with no co mparable activity for several 



tens of Myrs before and not much after (IKrabbe et al.lll995l ). Occam's razor would sug- 



gest that only one disk of WR/0 stars was present initially. The 10° sigma thickness of 



(Perets et al. 


2008 


), massive 


naries ( 


Perets et al. 


2008a 





of young stars at very high inclination s from the main clo ckwise disk structure, and/or the 



coherent counter-clockwise structure fICuadra et al.l l2008l ). Moreover, even more naassive 



and less likely perturbers such as in-spiraling intermediate mass black holes (lYu et all 120071 : 
Lockmann fc Baumgardtl l2008l ) have difficulties explaining both the clockwise and the in- 
clined counter-clockwise systems. In addition they require a fine tuned scenario, where the 
in-spiral occurs shortly after the stellar disk formation. It is therefore more likely that the 
clockwise and counter-clockwise configurations formed together and were not due to some 
dynamical evolution of a single disk-like structure. 

In situ star formation following t he infall of a single molecular cloud was recently studied 
in simulations (IBonnell fc Ricd 120081 ). In such a scenario a gaseous disk is formed and then 



fragments to form a single disk of stars (iNayakshin fc Cuadrall2005l : [Alexander et al.l 12008 



Bonnell fc Ricell2008l ). Efficient circ ularization in the original infalling gas lea ds to low orbital 
eccentricities of the massive stars ( ISanderd Il998l : IVoUmer fc Duschll 120011 ). A significant 
average eccentricity of the disk stars of 0.36 ± 0.06 pro bably requires that th e stars form on 
an orbital timescale, following the initial compression (IBonnell fc Ricell2008l ). However, the 
formation of warps in the disk and/or highly inclined stellar structures or even just a large 
fraction of stars outside the disk have not been observed in these simulations. 



Recent simulations by iHobbs fc NayakshinI (120081 ) suggest that strongly warped disks 
could form in-situ following cloud-cloud collisions in the GC. Moreover, highly inclined struc- 
tures such as stellar filaments, streams and small clusters could be observed in such simula- 
tions. Such scenarios could, in principle, naturally explain our observations, showing both 
a main warped disk and another highly inclined stellar structure, possibly with other less 
pronounced smaller structures and stars with various angular momenta outside the main 
structures. The second main structure could also be evidence of a second disk, possibly in a 
state of dissolution, or anot her coherent structure as seen in cloud-cloud collision simulations 
(IHobbs fc NayakshinI I2OO8I ). We caution, however, that such simulations, can currently ex- 
plore only a small fraction of the possible phase space for cloud collisions, and do not include 
radiative transfer effects. The present data set does not yet allow us to distinguish between 
a disk in a state of dissolution and stellar filaments, streams and small clusters 
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The eccentricities of the stars in the observed structures and the thickness of these 
structures could be explained both by the in-situ formation of ecceritric disks /structures 
jAlexander et aliboosl : ISonnell fc Ricelboosi : Iflobbs fc Navakshin|[2008l : IWardle fc Yusef-Zadeh 
20081 ) and/or through the dynamical heating of more circular configurations mainly by cusp 
stars and stellar mass black holes which excite the mean eccentricities to be as high as 0.4. 
The latter possibility, however, is not likely to excite eccentricities as high as 0.9, as possibly 
observed for some of the candidate disk stars (if they indeed belong to the disk). 

Another possibility that has been discussed in the literature is nuclear disk warping and 
eccentricity excitations, due to interaction with other massive non-spherical structures. Any 
configuration of two non-spherical symmetric populations causes the structures to precess. 
Differential precession at different speeds are observed as warps. There is a coherent rotation 
of the clust er of late-type stars in the Galactic Center but their distribution is spherically 
symmetric (ITrippe et al.l 120081 ) and hence does not cause any precession. Other planar 
structures in the Galactic Cente r such as the minispiral and the CND have too low masses 
( jJackson et al.lll993l : lLisztll2003l ) to cause significant warps. 



On the other hand, an in situ formation of the WR/0 stars in the Galactic Cen- 
ter in two disks, as possibly suggested by our observations, provides two massive inclined 
structures. Such a configuration has been suggested to induce w arping and eccentricity 



excitations through differential precessioii and a Kozai resonance (iNayakshin et al.l 12006 



Lockmann &: Baumgardtll2009l : ISubr et al.ll2008l ). The expected magnitud e of a warp due to 
differential precession induced by two disks can be estimated analytically. iNayakshinI (120051 ) 
assumes one massive ring (mass Mi) at radius Ri and a less massive test particle ring (mass 
M2) at radius R2. Let a be the angle between the angular momentum vectors of the two 
rings. The precession frequency cup of the less massive ring around the more massive one is 
than given by 



a;(i?2) _ 3 Ml 
n{R2) 



F?3 d2 



cos a 



4 M. {Rj + i?2)5/2 



GnM. 
3~ 

2 



^3-^ is the angular velocity of the less massive ring. For fixed masses the 



where fl = ^ = 

precession frequency has a maximum at -R2,max = v^3/7i?i. Therefore the precession time 
Tp is 



^ ^4M. 1 {Rl^.. + Rlf^\ jRl 

p > — — :r7^ 27r 



3 Ml COS a Ri^^^Ri 



GnM. 



With M, = 4 X 1^1 we get 



Tp > 10 years 



{Ri/?>"fl^ 
Mi/(5OOOM0) ■ 



(9) 



(10) 
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Paumard et al.l (120061 ) estimate the total mass of the counter-clockwise system to be 
SOOOMq. a counter-clockwise ring of about 5000 solar masses at a distance of about 3" from 
the Galactic Center, causes precession periods in the clockwise system of about 10'' years. 
This is compatible with a warp due to differential precession in the clockwise disk of about 
65°. This estimate is approximate. It is significantly modified for radially extended disks, 
or for stellar disks with eccentric orbits. A counter-clockwise ring with a radius of about 3" 
would cause a maximum precession of the clockwise disk at a radius of 3" and less precession 
for smaller and larger distances. In case of an initially flat clockwise disk, the innermost and 
the outermost border of that disk would have similar rotation axis, contrary to the data. This 
disagreemer it is intriguing but nee ds to be followed up with future more detailed models in 
the spirit of iNayakshin et al.l (120061 ) with updated constraints on the stellar orbits. Moreover, 
the above estimation of the warp amplitude did not take into account the influence of the 
cusp stars in which the disks are embedded. Analytical studies suggest that the stellar cusp 
around the disks would qu ench their strong i n teraction and would not allow for the Kozai 
resonance to be effective (llvanov et al.l l2005l : iKaras fc Subd 120071 : IChangj 120081 ). N-body 
simulations show that the cusp stars stabilize the disks and the eccentricity excitation is 
quenched (Perets et al., in prep.). If theoretical models with an initially fiat distribution 
cannot match the configuration seen in the data through warping, it would imply that the 
disks were born warped. This would provide an additional and very important constraint to 
the models. 

It is also interesting to note a well known example of another warped nuclear disk 
surrounding a supermassiye black hole in another ga laxy; the maser disk observation of 
NGC 4258 jMivoshi et al.lll995l : Ee^rnstein et al.1ll996h . The NGC 4258 disk extends from 
0.13 to 0.26 pc from the nucleus and shows a moderate warp of ~ 8° (iHerrnstein et al.lll996l ). 
Although interesting by itself, the observed nuclear maser disk shows much weaker warping 
than observed in the clockwise disk in the GC. It is likely to be of a dif ferent origin of the GC 
disk, and various suggestions have been discussed iri the literature (e.g. 



Arnaboldi fc Galletta 



19931 : ICaproni et al.ll2006l : iFerreira fc Ogilviell2008l ) ]Ulabav-Siddiki et all (120071 ) found that 
there are stable configurations of single massive nuclear disks in galaxies with pronounced 
warps for a wide mass range of the nuclear disk. 

The in-situ formation model requires that the gas densities in the infalling clouds are 
sufficiently high so that star formation occurs rapidly after the infall (but not before), so 
that the massive stars end up in the central parsec where they are observed and yet the gas 
does not have time to dissipate to one pla ne. Such a scenario produces a surface numbe r 
density of stars in the disk of S(r) oc (ILin fc Pringlelll987l : iHobbs fc NayakshinI l2008l ) . 
com patible with the obser vations. The top heavy initial MP expected from such a scenario 
(e.g. iBonnell &: Ricd 120081 ) is also consistent with our observations. We note, however, that 



- 49 - 



in some of the simulations by iHobbs &: NayakshinI (120081 ) they find that different regions in 
the main disk or in the other structures have a different mass function. The observation of 
more candidate member stars of the clockwise system might be required to confirm or refute 
this. 



6.2.2. Cluster Inf all 



In most scenarios for the infall of a cluster, the stars remain close to the in-spiral 
plane (IBerukoff fc Hansen! 120061 ) and form a si ngle flat star distribution with a thickness 
comparable to or larger than the original cluster. iLevin. Wu fc Thommed (120051 ) have shown 
that after the infall of a cluster of young stars around an intermediate mass black hole on 
an initially eccentric orbit a signiflcant fraction of stars get pushed out of the in-spiral 
plane and end up with large inclinations. However, the stars with large inclinations do not 
have a preferred orbital orientation; thus it is difficult for this scenario to form a second 
inclined coherent structure. Therefore, in the framework of the infalling cluster scenario, the 
observation of two disks of young stars requires two cluster infall events in the last 6 x 10^ 
years, rendering this scenario even more unlikely. 

In the infalling cluster scenario the stars which are stripped from the cluster as it 
spirals in have similar inclinations and eccentricities as the cluster itself requiring an initial 
eccentricity of the cluster of 0.36 ± 0.06. An in-spiral through dynamical friction would 
suggest a more moderate eccentricity. Nevertheless, obs ervations of the GC Arches cluster 
show that it is on a non-circular orbit (jStolte et al.ll2008l ): although at a large distance from 
the Sgr A*. 

The observed warps give strong evidence against the infalling cluster scenario. In the 
framework of the infalling cluster scenario, the stars will keep the mean angular momentum 
orientation of the initial cluster if the cluster sinks in rapidly resulting in a fiat, non- warped 
disk of stars. To form a warp of the observed magnitude through differential precession, the 
stars need to have been in a relatively thin disk for a long time. But in the infalling cluster 
scenario, the stars barely had the time to get to their sub-pc orbits, and thus they could not 
have spent much time disk. 

The infal ling cluster scenario predicts surface density profiles a s shallow as S(r) oc 

^-0.75 j^ggg g_g_ 

iHansen fc Milosavljevi^booal: iBerukoff fc Hansen! 120061 ). This is considerably 



shallower than that observed. 



Giirkan fc J2005h argue that this discrepancy 



can, m 



principle, be overcome by initial mass segregation in the cluster. The lower mass stars of 
the cluster are lost at larger distances and the most massive stars inside. This proposition 



will be tested by a fu ture search for lower mass B stars at distances beyond 12" from Sgr A* 
flMartins et aDboOoi ). 



6.3. Summary 

To sum up, our observations are compatible with an in situ formation of these WR/0 
stars in a clockwise rotating disk and another highly inclined counter-clockwise structure, 
possibly a disk. Although the infall of two clusters may possibly explain the two systems of 
WR/0 stars, the stars would have remained in a disk for too short to develop the observed 
warp. 

7. Conclusions and Outlook 

Our latest sample of early-type stars between 0.8" and 12" from the Galactic Center 
contains 90 WR/0 stars. The most important results are: 

• 55% of the WR/0 stars are candidate members of a clockwise disk 

• the clockwise disk shows a large warp (in total (64 ±6)°, at a significance level of > 10 
sigma), the disk angular momentum direction is a function of the projected distance. 

• the clockwise disk is not locally fiat, the inclinations of the clockwise system are com- 
patible with a two-dimensional Gaussian distribution with a sigma of 10° 

• the clockwise system is not circular; it has a mean eccentricity of 0.36 ± 0.06. 

• 20% of the WR/0 stars are candidate members of a coherent feature amongst the 
counter-clockwise stars (pre-trial significance of 4.5 sigma for the interval of projected 
distances between 3.5" and 7", corresponding to a 2% chance that this feature is a 
statistical fluctuation in an isotropic cusp star distribution). 

• the counter-clockwise system may be a disk in a dissolving state: 10 out of 29 counter- 
clockwise rotating stars have angular separations larger than 30° from the counter- 
clockwise system, while only 11 out of 61 clockwise moving stars have a separation 
larger than 30° from the clockwise disk. 

The observation of a warped disk of WR/0 stars in the Galactic Center has substan- 
tially strengthened our previous conclusion that the population of 100 young, massive stars 
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arranged mainly in two coherent structures between 0.8" and 12" from Sgr A* is due to in 
situ star formation (as opposed to an infalling cluster scenario) from gas that fell into the 
central region about 6 x 10^ years ago and became dense enough to overcome the tidal forces. 



A. Coordinate systems 



As in ( iPaumard et al.ll2006l ) we used the usual Cartesian coordinate system in offsets 
from Sgr A*: x = cos5da increased eastward, y = d6 increased northward and z = dD 
increased along the line of sight away from the observer. We used a spherical coordinate 
system (x, y, z) (r, 9, 0): 



X = r sin 6 cos i 
y = r sin 6 sin <; 
z = rcosO . 



(Al) 
(A2) 
(A3) 



Figure [25] illustrates the definition of the classical orbital elements used in this work in 
the framework of the chosen coordinate system in this work. 
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